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Comparison  of  Two  1550  nm  Ultra  Narrow-Band  Optical  Infinite 
Impulse  Response  Filters  for  High-Speed  Optical  Signal  Processing 

Shyh-Lin  Tsao*  and  Te-Yu  Chen 
Department  of  Electrical  Engineering,  Yuan  Ze  University 
Chung-Li, Taiwan,  R.O.C. 

ABSTRACT 

Comparison  of  the  frequency  responses  of  two  high-speed  optical  fiber  processing  elements  is  reported. 
The  two  filters  we  considered  are  a two-loops  optical  infinite  impulse  response  (HR)  filter  and  a three-loops 
optical  HR  filter.  The  theoretical  models  of  these  two  kinds  of  optical  HR  band  pass  filters  are  built  up. 
Some  suitable  coupling  coefficients  are  searched  for  simulating  the  narrow  band-pass  filters.  The  frequency 
responses  under  different  conditions  are  analyzed.  The  theoretical  models  can  provide  parameter  evaluation 
for  experimental  design.  A higher  design  flexibility  can  be  given  by  the  three-loops  optical  HR  filter.  These 
filters  can  be  implemented  with  PIC  (photonic  IC)  techniques  and  applied  for  optical  data  signal  processing. 

Keywords  • Optical  Data  Storage,  Optical  Signal  Processing,  Optical  IIR  filter,  Optical  Fiber 
Communication 


l.INTRODUCTION 

In  the  next  generation  memory  system,  high  capacity  optical  data  storage  will  be  the  future  trend.  Many  methods,  such 
as  hologram,  CD-R,  DVD,  are  studied  for  high  capacity  data  storage  above  Giga-byte 1-3  . In  recent  years,  some  fast 
data  access  and  optical  signal  processing  methods  are  extensively  reported4  8 . Especially  reading  optically  storaged 
data  by  fiber,  then  using  time-space-conversion  optical  signal  processing  attract  a lot  of  attentions 5 . 

Promising  progress  is  shown  in  high  speed  optical  communication  systems  developments  recently 9 ,8  . The  fiber  optic 
recirculating  delay  line  is  one  of  the  basic  optical  signal  elements  generally  being  studied.  We  also  reported  some 
various  designs  and  applications  for  high  speed  signal  processing19-21  . Ultra  fast  signal  processing  over  100  Gbps  is 
difficult  to  do  by  traditional  electric  signal  processing  circuit.  This  optical  digital  signal  processor  (OPSP)  may  be 
applied  in  high-speed  DVD  data  storage  systems  as  shown  in  Fig.l.  In  this  paper,  we  design  and  analyze  two  ultra 
narrow-band  optical  infinite  impulse  response  IIR  filters  for  high-speed  optical  signal  processing.  Those  two  designed 
narrow-band  optical  IIR  filters  are  operating  at  wavelength  of  1550nm.  The  center  of  the  narrowband  filter  is  located 
at  1 00GHz.  We  use  Z-trausform  to  analyze  the  frequency  responses  of  these  two  IIR  filters : two-loops  optical  IIR  filter 
and  three-loops  optical  IIR  filter.  The  two  theoretical  models  for  comparing  their  frequency  responses  are  built  up. 
Optimal  coupling  coefficients  are  sorted  for  simulating  narrow  bandwidth  optical  band  pass  filters. 
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2.  THEORETICAL  MODELS 

We  use  Z-transform  techniques  to  build  up  the  theoretical  model  of  a two-loops  optical  HR  filter  and  a three-loops 
optical  HR  filter,  respectively  in  this  section.  Assume  these  two  kinds  filters  are  linearly  time-invariant  systems  and  no 
temperature  fluctuation  effect,  we  derive  the  transfer  functions  as  following. 


2.1  Two-Loops  Optical  HR  Filter 

The  schematic  diagram  of  two-loops  optical  HR  filter  with  two  3*3  optical  couplers  and  four  optical  delay  lines 
is  shown  in  Fig. 2.  In  Fig.3,  we  modeled  the  equivalent  theoretical  model  by  signal  flow  diagram.  The  signal  can  be 


represented  as  Eq.  (1)  to  Eq.  (6) 

x,  (Z)  = K3  [I,  (Z)  + rX4  Z'XtXA(Z)  + rXj  Z'Xj  X3(Z)]  (1) 

X2  (Z)  = K2  [I,  (Z)  + rX4  Z'XjX4(Z)  + rX3  Z'x>  X3  (Z)]  (2) 

X3(Z)  = K6[I2(Z)  + rX2  Z'X2X2(Z)  + rXi  Z‘x'  X,  (Z)]  (3) 

X4(Z)  = K5[I2(Z)  + rXjZ'XjX2(Z)  + rXiZ‘XlX[  (Z)]  (4) 

o,  (Z)  = Kj  [I,  (Z)  + rX4  Z-Xj  X4  (Z)  + rXj  Z'x>  X3  (Z)]  (5) 

O 2 (Z)  = K 4 [1 2 (Z)  + rx  2 Z'Xj  X 2 (Z)  + rXj  Z‘x‘  X j (Z)]  (6) 


I,  (Z)  and  I2(Z)  are  the  input  signals  at  the  input  ports.  0,(Z)  and  02(Z)  are  the  output  signals  at  two  ports.  X,(Z), 

X2(Z),  X3(Z),  X4(Z)  represent  the  input  signals  of  delays  Z'X| , Tx- , Z'Xj , Z'x> , respectively.  The  coupling 
coefficients  and  path  transmittive  coefficients  for  each  optical  output  are  represented  as  Kj(j=1_6)  and  rxj(j=1_6)=l, 
respectively. 

The  relationship  of  the  above  six  equations  can  be  simplified  as 

X,  (Z)  = K3[I,  (Z)  + ^ rXj  Z‘X4X3(Z)  + rX3Z‘x’  X3  (Z)]  (7) 

X3  (Z)  = K3[I2  (Z)  + ^2-  rX2  zx2  x,  (Z)  + rXlZ‘x'  X,  (Z)]  (8) 

K3 

Substitute  Eq.  (8)  into  Eq.  (7),  we  find 

X|(z)  K3[li(z)t(K,rX4Z-x«  tKt,x,Z-xih(z)]  (?) 

and 


A = 1 ■ 


K3K6: 


X,rX3 


z-(x,+x3)  _K3K5rXIrX4Z-(x‘+X4>  -K2K6rXirX3Z-'X2+X2)  -K2K5rX2rX4Z 


-(X2+X4) 


Similarly,  we  can  write  X3  (z)  as 


(10) 
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(11) 


, 1 K* 
V I7I_  ......  _ 

I2(z)+I 

(K3rxiZ-x>  +K2rX2Z-x= 

Mz) 

A 

Then  the  output  signals  can  be  represented  as 
„ K,[l,(z)+(K5rX4Z-X4  + K6rXjZ_X3 ) 

i2(z)] 

0 (Z)-^ 

[l2(z)+( 

A 

K3rxiZ-x-+K2rx?Z-x=) 

>i(z)l 

A 

(12) 

(13) 


In  this  paper,  we  assume  no  input  signal  launched  to  port  I2.  Therefore,  the  above  two  transfer  funtions  can  be 
simplified  as 


H„(Z)  = 


9M 

h{z) 


Ki 

A 


H ,2(Z)  = 


Q](z)  K4(K3ry2-x'<-K,r>iZ-’<-) 
I,(Z)  A 


(14) 


(15) 


2.2  Three-Loops  Optical  HR  Filter 

The  schematic  diagram  of  a three-loops  optical  HR  filter  is  shown  in  Fig.4.  This  filter  includes  two  4x4  optical 
couplers  and  six  optical  delay  lines.  The  schematic  diagram  of  theoretical  model  of  this  filter  is  illustrated  in  Fig.5. 
Similarly,  we  consider  the  nodal  signals  X,(Z),  X2(Z),  X3(Z),  X4(Z),  X5(Z),  and  X6(Z),  represent  the  input  signals  of 

the  phase  delays  Z'X| , Z'x- , Z‘Xj , Z"X4 , Z"X5 , Z'x° , respectively.  The  transmittive  coefficient  of  each  delay 
path  rXj(j=,_8)is  set  to  be  1.  The  coupling  coefficient  of  each  optical  output  is  represented  as  which  is  a positive 

real  number.  Because  of  energy  conservation,  the  signal  in  each  coupler  should  satisfy  K^Kj+Kj+^+Ks+K^Ky+Kg 
= (1-L),  and  L means  the  intrinsic  loss  of  optical  couplers.  In  this  paper,  we  assume  the  intrinsic  loss  is  small  enough 
that  can  be  neglected.  According  to  Fig.5,  the  transfer  functions  can  be  derived  according  to  the  following  equations  as 


X , (Z)  = K 4 [I , (Z)  + rx  ^ Zx<  X 4 (Z)  + rXs  Z'A>  X 5 (Z)  + rX6  Z‘*‘  X 6 (Z)] 


-x. 


7-X, 


7‘X, 


X3(Z)  = Kj  [I,  (Z)  + rX4  Z‘Xj  X4  (Z)  + rX;  Z’x’  X5  (Z)  + rX6  Z'x‘  X6  (Z)] 
X 4 (Z)  = Kg  tl  2 (Z)  + rX|  Z'X|  X , (Z)  + rXj  Z'x>  X 2 (Z)  + rXj  Z'x’  X 3 (Z)] 
X 5 (Z)  = K 7 [1 2 (Z)  + rx  Z'x'  X , (Z)  + rx  X 2 (Z)  + rXi  Z x'  X 3 (Z)] 


X6 (Z)  = K6 [I2 (Z)  + rx  Z’A'  X,  (Z)  + rx  Zx> X2  (Z)  + rx,  Z^X3 (Z)] 


Z'X 


,-x4- 


7*X< 


z-X, 


6X6(Z)] 

(16) 

;‘x6(Z)] 

(17) 

:‘X6(Z)] 

(18) 

5X3(Z)] 

(19) 

JX3(Z)] 

(20) 

JX3(Z)] 

(21) 

x6(Z)] 

(22) 
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(23) 


02  (Z)  = K5[I2(Z)  + rX|  Z'x'  X,  (Z)  + rX2  Z‘x>  X2  (Z)  + rxZ^X3(Z)l 

Similarity,  we  deduce  the  above  equations  and  find  that 

X,  (Z)  = K4[Ij  (Z)  + rX4Z-x<X4(Z)  + rx$  -~7-  Zx<  X5(Z)  + rX6  ^ Z'x'X6(Z)] 

X 4 (Z)  = K 8 [1 2 (Z)  + rx , Z'x'  X , (Z)  + rXj  ^ Z'x‘  X 2 (Z)  + rx  3 ^ Zx'  X 3 (Z)] 

K.4 

From  the  above  two  equations,  X,(Z)  can  be  derived  as 

K 4 [l,(z)+(K8rX4Z-x-  + K7rX;Z~X;  + K.r^Z^  )l2(z)] 

A' 


X,(Z)  = - 


and 


A'=l-K4K8rXirX4Z-(x'+x^ 

7-(X2+X5) 


K4K7rXlrx  Z-(X‘+X^  -K4K6rx  rx  Z^5™  - K3K8rx  rx  Z^^ 


■K3K7rX2rX5Z-^+A5)  -K3K6rX2rX6Z-(5™  -K2K8rXirX4Z-(X’+x^  -K2K7rX3rX5Z-(X’+x’> 

rx./^ 


•K2K6rX5rX6Z-(X’+x‘) 


The  X4(Z)  can  also  be  derived  as 

K8[l2(z)  + (K4rx  Z"Xl  +K3rXjZ-x>  + K2rXjZ-x’ 


X4(Z)  = - 


A’ 


After  a lot  of  reductions , we  can  find  the  transfer  functions  of  the  three-loops  HR  optical  filter  as 

Ii(z)  A' 


H12(Z)= 


02  (Z)  K5  (K4rX)  Z_x'  + K3rX2  Z‘x>  + K2rXlZ~x» ) 


i,(z) 


A' 


(24) 

(25) 

(26) 
(27) 

(28) 

(29) 

(30) 


3.  SIMULATION  RESULTS  AND  COMPARISONS 

Here  we  use  the  theoretical  model  derived  in  the  above  section,  the  frequency  responses  of  these  two  optical  HR  filters 
are  simulated.  The  simulation  results  can  be  applied  for  practical  realizing  various  filters.  Taken  narrow  bandwidth 
optical  band  pass  filters  as  example,  the  optimal  parameters  searching  work  is  performed. 

Table  1 shows  three  set  of  searched  parameters  corresponding  to  the  simulated  three  examples  of  frequency  responses 
of  the  two-loops  optical  IIR  filter.  According  to  the  theoretical  model  we  derived  in  subsection  2.1,  three  examples  of 
frequency  responses  of  two-loops  HR  filters  are  shown  in  Fig.6.  Example  I is  the  best  result  we  found.  If  the  parameter 
K4  increases,  the  sidelobe  will  arise  as  shown  in  the  curveture  of  example  II.  The  delay  coefficients  X,,  X2,  X3,  X4  can 
be  adjusted  for  designing  higher  order  filters.  The  simulation  result  example  III  shown  in  Fig.6  give  a general  scope  of 
view  of  higher  order  filter. 

The  parameters  we  used  for  simulating  three-loops  IIR  filters  in  Fig.7  are  listed  in  Table  2.  Similarily,  the  numerical 
results  are  shown  as  example  I,  II  and  III  for  comparing  the  results  of  Fig.6.  The  optimal  result  of  narrowband  filter  is 
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shown  in  example  I.  When  K,  increases,  the  sidelobe  will  grow  up  as  shown  in  example  II  like  the  two  loops  HR 
filters.  Increasing  the  order  of  this  filter,  the  ripple  of  this  filter  will  increased  as  shown  in  example  III.  We  found  the 
ripple  is  larger  then  the  two  loops  HR  filters  if  the  parameter  is  not  well  choosen. 

Comparison  of  the  optical  two  numerical  results  in  Fig. 6 and  Fig.7,  we  zoom  in  the  two  cases  of  example  I and  show 
the  spectra  details  in  Fig. 8.  We  can  find  that  the  three-loops  HR  filter  has  a narrower  optical  frequency  bandwidth  than 
the  two-loops  HR  filter.  The  narrowest  bandwidths  of  the  two-loops  IIR  filter  and  the  three-loops  IIR  filter  we  found 
are  0.00012  Hz  and  0.00007  Hz,  respectively.  The  narrow  linewidth  can  provide  pure  filtering  effect  for  extracting  the 
phase  lock  loop  carrier  signal  in  DVD  control  systems. 

4.  CONCLUSIONS 

The  access  data  rate  of  optical  storage  systems  increases  rapidly  in  the  last  five  years.  Therefore,  the  signal  processing 
will  follow  the  trend,  increasing  the  processing  speed  in  the  future.  Accompany  with  the  age  of  optical  data  storage 
above  Tera  byte,  studying  all-optical  fiber  high  speed  signal  processing  elements  is  important.  The  research  results  in 
this  paper  are  very  helpful  for  practical  realizing  various  optical  signal  processing  filters.  The  narrow  band  filter  at 
100GHz  can  be  applied  for  clock  recovery  in  high  speed  data-accessed  DVD  systems  accompany  some  fiber  pick  up 
heads. 
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Table  1 


K, 

k2 

k3 

k4 

k5 

k6 

X, 

X2 

x3 

x4 

I 

0.0001 

0.49995 

0.49995 

0.001 

0.4995 

0.4995 

1 

2 

1 

2 

n 

2 

4 

4 

2 

4 

4 

1 

2 

1 

2 

0.0001 

0.49995 

0.49995 

0.9 

0.05 

0.05 

3 

4 

5 

6 

